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ABSTRACT 



The formation of CO2 in quiescent regions of molecular clouds is not yet 
fully understood, despite C0 2 having an abundance of around 10-34 % H 2 0. 
We present a study of the formation of CO2 via the non-energetic route CO + 
OH on non-porous H2O and amorphous silicate surfaces. Our results are in the 
form of temperature-programmed desorption spectra of CO2 produced via two 
experimental routes: 2 + CO + H and O3 + CO + H. The maximum yield 
of CO2 is around 8 % with respect to the starting quantity of CO, suggesting 
a barrier to CO + OH. The rate of reaction, based on modelling results, is 
24 times slower than O2 + H. Our model suggests that competition between 
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CO2 formation via CO + OH and other surface reactions of OH is a key factor 
in the low yields of CO2 obtained experimentally, with relative reaction rates 
kco+H < kco+OH < k H2 o 2+H < k H+H, ko 2 +H- Astrophysically, the presence of 
C0 2 in low Ay regions of molecular clouds could be explained by the reaction 
CO + OH occurring concurrently with the formation of H2O via the route OH 
+ H. 

Subject headings: astrochemistry — ISM: molecules — methods: laboratory 



Introduction 



The first observations of solid CO2 (henceforth CO^Q were made by the InfraRed 



Astronomical Satellite (IRAS. iD'Hendecourt fc Jourdain de Muizon I (119891 ) ) . The molecule 



(de Graauw et al. 


1996 


), massive protostars ( 


Gerakines et a 


. 1999: 


Gibb et al. 


2004 


), low 


mass YSOs ( 


Nummelin et al. 


2001; 


PontODDidan et a 


. 12008 


), background stars ( 


Knez et al. 


2005). and in other galaxies ( 


Shimonishi et al. 


2010|; 


Oliveira et al. 


2011). Based on these 



observations, C0 2 ( s ) is seen to be seemingly ubiquitous, and one of the most abundant solid 
phase molecular species, approxim ately 10-34 % H 2 Q. It is be lieved to form in the solid phase, 



due to low gas phase abundances (Ivan Dishoeck et al. 



1 996 ), with evidence suggesting that 



much C0 2 ( s ) production occurs in quiescent regions (IPontoppidan 
20011 ). yet the key question remains: how does C0 2 ( s ) form? 



2006; Nummelin et al. 



Many experimental studies have been perform ed to study the energet ic formation routes 
to CQ 2 . Irradiation of pure CO ices wit h photons (jGerakines et al. II1996I) . charged particles 



( jPalumbo et al. 1 119981 ) , and electrons ( jJamieson et al. 1 120061 ) h ave yielded CO2. Similar 



experiments with mixtures of CO and H 2 Q were also successfu l (jEhrenfreund et al. 111997 



Palumbo et al. 1 11998c lloppolo et al. I l2009t lLaffon et al. 1 120101 ) . The irradiation of hydro 



genated carbon grains with ions and electrons produced small quantities of CO and C0 2 
rtMennella et~al1l2004i liooS i. 



C02( s ) is abundant in quiescent, as well as star-forming re gions. While the ro le of 
energetic pathways can not be discounted entirely in these regions ( jWhittet et al.lll9981 ). the 
study of non-energetic formation routes is fundamental to fully understanding the obse rved 
abundances of C02( s ). Potential non-energetic formation routes (jRuffle fc Herbstll200ll ) are: 



CO + O -> CO, 



(1) 



CO + OH -> C0 2 + H 



(2) 
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ECO + O ->■ C0 2 + H 
HCO + 0H ->■ C0 2 + # 2 



(3) 
(4) 



Theoretical studies of route (jT| suggest that the formation of C0 2 proceeds with a 
high barrier of around 2500 - 3000 K, lowered on surfa ces via the hot O atom or Eley-Rideal 
mechanisms (ITalbi et al. 1120061 ; iGoumans et al. II2008I ). A solid phase study determined that 
this pat hway was feasible o nly via reaction in water pores, under a water ice cap, and upon 
heating (IRoser et al. 1120011 ). suggesting that it would not occur under the conditions present 
in quiescent molecular clouds. Reactions fl3]) and (TjJ have never been studied expressly in 
the solid phase. 



Route (El) has been extensively studied in the gas phase, both e xperimentally (|Frost. Sharkey fc Smith 



1993 



2001 



Fulle et al. 



1996 



Chen fc Marcus 



2005 



Baulch et al. 1 120051) a nd theoretically (lYu. Muckerman fc Sears 



Sun &: Law 1120081 ). due to its importance in atmospheric and 
combustion chemistry. Recently, reaction ([2]) was experimentally studied in the solid phase 
by Refl ection Absorption InfraRed Spectroscopy (RAIRS) for the first time, with positive 
results (jOba et al. Il2010l ). Due to the limitations of the adopted method, it was not possible 
to produce a pure beam of OH, and therefore the chemistry is difficult to constrain with 
a simple series of reactions; in particular, it was experimentally complex to distinguish be- 
tween (J2J) — (SI) - OH was produced in the gas phase via a plasma discharge of H 2 0, a process 
which yields a mixture of products including OH, H, H 2 , O and 2 . Although it is claimed 
that all OH radicals are in the rovibrational ground state due to collisions with the beam 
walls, well defined spectroscopic studies of plasma discharges suggest that interaction with 
the walls is likely to lead to OH recombination, rather than yield ground state OH, and 
that th e major component s of a plasma of H 2 Q are H 2 and H 2 0, with lower abundances 
of OH (IMedard et al.ll2002l ; iFujii et al. II2002I ). Furthermore, experiments in the absence of 
CO produced H 2 2 and 3 , whose yields varied with surface temperature, suggesting that 
surface temperature itself, mobility of the discharge products on the surface, and, poten- 
tially, the desorption rate of CO from the surface, rather than the rovibrational state of 
the OH, is responsible for the changing yields of C0 2 observed at different temperatures. 
Finally, reaction fl2]) was found to proceed with little or no barrier, suggesting the presence 
of rovibrationally excited OH. We contend that, under these conditions, the C0 2 yield can 
not be assumed to be independent of the excitation state of OH and thus further study of 
reaction (j2J) is imperative. A subsequent RAI RS study produced O H in the solid phase from 
a mixture of 2 :CO in a multilayer regime (lloppolo et al.l 1201 laf ). Due to the multilayer 
regimes investigated, both previous studies also involved more complex chemistry than sim- 
ply C0 2 formation, such as the formation of H 2 C03 and CH 3 OH, which further complicate 
the quantitative analysis of reaction (j2J). 
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Here we present the first temperature-programmed desorption (TPD) spectra of C02( s ) 
formed via ([2]) in the solid phase, under interstellar conditions of temperature and pressure. 
In this study, OH was produced on the surface by reaction of 2 and 3 with H in order to 
constrain better the reaction pathways in the system. The reaction was studied on both an 
amorphous silicate and a non-porous water surface, in a low coverage regime, with the aim of 
limiting chemistry to only CO2 production. In contrast to previous studies, a simple kinetic 
model was developed to determine relative reaction efficiencies and calculate the activation 
energy of reaction (T2J). 



2. Experimental 



Experiments were performed using the FORMOLISM apparatus (lAmiaud et al. II2006I ). 
Briefly, the experimental set-up consists of an ultra high vacuum (UHV) chamber (base 
pressure ~ 10 ~ 10 mb ar), containing an amorphous silicate-coated copper surface (5-400 K, 



Lemaire et al.l ( 120101 )). Molecules are dosed onto the surface via two triply differentially 
pumped beam lines. Desorption of molecules from the surface is monitored using a quadrupole 
mass spectrometer (QMS), positioned directly in front of the surface. Experiments were 
performed on either bare silicate, or a non-porous, amorphous water film (nj?-H 2 0) of 
~100 monolayers (ML) grown on the silicate by spraying water vapour from a microchannel 
array doser (held at 120 K during water desorption, then cooled to 10 K before commencing 
the experiments). 

Two different surfaces were investigated in order to, firstly, mimic two interstellar en- 
vironments and, secondly, to determine the surface dependency of route (j2D- Amorphous 
silicate is an approp riate mimic of int erstellar dust grains, composed of siliceous and car- 
bonaceous material (IGreenberg 1 12002| ). In molecular clouds these grains are covered in 
an ice mantle, the larges t component of which is H2O, at abundances of up to 100 ML 
( jWilliams fe Herbstl 120021 ). np-H 2 was used in this study to eliminate the complexity of 
chemistry occurring in pores. 

Neither O3 nor O2 has been observed in an interstellar ice, so these experimental condi- 
tions are not directly astrophysically relevant, but were used to produce OH in a controlled, 
reproducible manner. It is experimentally complex to create, maintain, and deposit onto a 
surface, a pure, stable beam of the OH radical in the ground state. Thus, in this work, OH 
was produced on the surface via two routes: the hydrogenation of O2 and of O3. O2 is easier 
to utilise experimentally, but a study involving both species constrains reaction mechanisms 
better than using a single species. Due to the limits of sensitivity of the QMS, quantities 
below 0.1 ML were not investigated. 



-5 - 



OH was not measured directly on the surface, due to its short lifetime, but the produc- 
tion of O2, H2O2 and H 2 0, during control experiments on the hydrogenation of O2 and O3, 
confirms its presence, according to the reaction scheme: 

3 + H ^0 2 + OH (5) 
2 ^ H0 2 ^ H 2 2 ^ H 2 + OH (6) 

Figure [1] describes the chemical network of OH relevant to these experiments. The OH 
radical, produced by (jSJ) and (jSJ), could react with CO as in (j2J) to produce CO2, or with H 
to form H 2 0, thus, the relative rate of these reactions is an important factor determining the 
yield of CO2 and the aim of this work was to elucidate the relative rates of these reactions. 

All experiments are summarised in Table [T| approximately 1.5 ML of O3 or 0.5 ML of 
O2 were dosed onto the surface via one molecular beam, followed by ~0.5 ML of isotopically 
labelled 13 CO via a second beam. Finally, H atoms were deposited, via a plasma discharge 
of H 2 on beam 1, for a range of exposure times of between and 20 minutes. The deposition 
rate of H on the surface was ~ 5 x 10 1 2 atoms cm~ 2 s~ 1 , tak ing into account the dosing 



pressure and the dissociation rate of H 2 (lAmiaud et al. 1 120071 ). During O3 deposition, the 
surface was held at 45 K to ensure that any traces of O2 p resent in the O3 beam desorbed 



from the surface (for detailed 3 production method, see iMokrane et al. I (120091 ) ) ; it was 
then cooled to 10 K before continuing. For all other molecules, the surface was held at 10 K 
during dosing. 

H atoms were hot when produced in the plasma discharge of H 2 , but cooled to room 
temperature before exiting the molecular beam, due to collisions with the walls. Isotopically 
labelled 13 CO was used to avoid contamination of the results by 12 C02 or 12 CO, pollutants 
present at very low gaseous concentrations in the chamber. 

To measure the products of the reaction, the surface was heated from 10 to 100 K; 
desorbing molecular species were monitored with the QMS. Each TPD cycle lasted approx- 
imately three hours. 



3. Results and Discussion 

Figure [2] shows TPD spectra of 13 C02 produced by H irradiation of O3 with 13 CO, 
and O2 with 13 CO. 13 C02 was produced during all experiments where H irradiation was 
performed. The 13 C02 desorbs from each surface over a similar temperature range, but 
with a slightly different peak shape, indicative of the roughness of the underlying surfaces. 
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It is clear that, since all data in Figure |2] are measured at the same H irradiation time, 
substantially more 13 CC>2 is produced by the reaction of O3 than O2. This is evident from 
reactions (JSJ and three hydrogenation steps are required to generate a single OH radical 
from O2, whereas O3 also generates an OH directly. 

As the solid blue line in Figure |2] shows, if the experiment was conducted without 
H irradiation, no 13 C02 was produced. Nor was 13 C02 production seen during control 
TPDs of 13 CO, O2 or O3. When 13 CO was not present, H 2 and H2O2 formed, rather 
than 13 C02- It was assume d that all the 13 CQ2 formed during H ir radiation, in agreement 
with previous experiments (lOba et al. I l2010t lloppolo et al.l l2011al ). Within the limits of 
measurement, no gas phase 13 C02 was observed during the irradiation, indicating that the 
13 C02 remained on the surface. No additional 13 C02 was produced during the TPDs because 
control experiments, where 3 was irradiated with H before deposition of 13 CO, yielded 
negligible 13 C02, suggesting that when OH is produced, it reacts quickly on the surface. 

When searched for, we saw no evidence of the production of H 13 COOH, validating 
the hypothesis that investigating a submonolayer coy e rage restricts the ch emistry to CO2 
production, unlike previous studies ( lOba et al. Il2010t lloppolo et al.l l2011ai ). Neither were 
H^CO nor 13 CH30H seen in desorption. However, H2, 12 CO, 13 CO, 12 C02 were seen during 
all TPDs; 13 C02 was seen upon H irradiation; H2O2, H2O were seen upon H irradiation, 
during extended TPDs to higher temperature; O3 was seen only when O3 had been deposited; 
O2 was seen in all experiments, except those with O3 and no H irradiation. 

Figure [3] illustrates 13 C02 production as a function of H irradiation time (solid symbols). 
It is clear that 13 C0 2 was produced at comparable rates on both surfaces (solid triangles 
versus solid squares), regardless of the starting material (O2 or O3), vindicating our earlier 
conclusion that the underlying surface does not play a significant role in this reaction. In 
every experiment, some 13 CO desorbed during the TPD, suggesting it was present in excess, 
never completely reacting with OH or H. In addition, as O3, O2 and H were dosed via one 
beam, while 13 CO was dosed via a second, even after alignment the maximum overlap at- 
tainable is less than 100 %, so not all reagents were dosed on the same region of the surface. 
Conversely, as Figure [3] shows, the reactions are very sensitive to the quantity of OH gener- 
ated, which itself depends upon the starting quantity of 2 or 3 on the surface. Although 
surface coverages were controlled to within ± 0.2 ML between experiments, this difference 
was sufficient to account for the varying concentration of 13 C02 observed in Figure [3J 

The maximum yield of 13 C02 was ~ 8 % with respect to 13 CO (Table [TJ experiments 
F,K); the presence of a complex barrier helps to explai n this. Gas phase and t h eoret i- 



cal studies predict a three-stag e barrier (~ 500 K, e.g. iFrost. Sharkey &: Smith I ( 119931 ); 



Yu. Muckerman Sears I (120011 ) ) . The reaction proceeds via an energetic HOCO interme- 
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diate, which isomerises from trans-HOCO t o cis-HOCO, befor e dissociating t o form CO2 



flSmith fc Zellner Ill973l : lAlada et al. Ill993l : iLester et al. Il2000h . iLester et al. I (120011 ) sug- 



gest that a precursor OH-CO complex forms prior to the HOCO intermediate. On a surface, 
the reaction probability is even more reliant upon the relative orientation of CO and OH. A 
recent theoretical study on a corone ne surface shows that O H physisorbs with the hydrogen 



atom pointing towards the surface (IGoumans et al. II2008I ). Compared to the gas phase 



the activation barrier to trcms-HOCO formation is slightly lowered, and the intermediate 
stabilised. A barrierless reaction between this stabilised HOCO complex with a further H 
atom could produce CO2 + H 2 . However, exp eriments suggest this reaction could also yield 
HCOOH or H?Q + CO Jloppolo et al.l boilbh . while reaction with OH could yield H 2 C0 3 



( lOba et al. II2010I ). From the results presented here, it is not possible to determine whether 



13 C02 formed from cis-HOCO — > CO2 + H, or by hydrogenation of HOCO, although no 
H 13 COOH was observed when TPDs were run to 200 K, thus we assume that, under our 
experimental conditions, reaction fl2]) produces only 13 C02- 

This complex barrier somewhat explains the low 13 C02 yields, but there are further 
constraints to be considered. Due to the low coverages investigated here, the probability of 
13 CO and OH meeting on the surface is small. OH recombination co uld produce H 2 2l o r 
the competitive reaction OH + H could remove OH from the surface ( jloppolo et al. II2008I ). 
Also, in these experiments, 13 CO was dosed after O3 or O2, so at high enough coverages 
it could block H from reaching these reagents, allowing C O hydrogenation to ar tificially 



dominate. Previous s tudies show that h ydrogenation of O3 (IMokrane et al. 1120091 ) and O 



( jloppolo et al. 



proceeds via: 



2008t iDulieu et al.ll2010l ) occurs with no barrier, while hydrogenation of CO 

(7) 



CO HCO ^ H 2 CO ^ H 3 CO CH 3 OH 



with a barrier of 390 K at 12 K to CO + H (IFuchs et al. M2009I : IWatanabe fc Kouchil 120021 ). 
Here, the 13 CO surface coverage was always below 1 ML, and neither H^CO nor 13 CH30H 
desorbed during extended TPDs, indicating that little 13 CO hydrogenation occurred. If 
H 13 CO was not produced in significant concentrations via ([7]), it follows that 13 C02 was 
not produced at measurable quantities via (jl]), not least due to the low probability of any 
H 13 CO produced encountering OH on the surface. These conclusions indicate that in this 
experiment 13 C02 formation occurred exclusively via reaction ([2]). 

A simple kinetic model of the experimental system was developed to describe the pro- 
duction of 13 C0 2 via ([2]), based on a set of coupled first order rate equations. Figure [TJ shows 
the potential detailed reaction scheme for these experiments. However, as illustrated above, 
given the low surface coverages employed here, a number of reactions, for example OH re- 
combination, can be eliminate d and are thus shown in grey in the figure. Neither CO2 nor 
HCOOH react further with H ( IBisschop et al.l 120071 ). so such reactions were also ignored in 
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the model. As discussed above, H 13 CO was likely not produced at significant concentrations 
in these experiments. However, the rate of CO hydrogenation under present conditions was 
constrained by control experiments with only CO on the surface. It was assumed that the 
underlying surface had no effect on the reaction rate (as illustrated by Figure [3]), so the 
model treats both surfaces simultaneously. The best fit to the experimental data was found 
by varying the rates of reaction, kj, of CO + OH and H2O2 + H, while constraining all other 
kj with empirical data (reactions shown in black in Figured]). That only two free parameters 
{kco+OH and k# 2 o 2 +ii) were required to fit all of the data presented here, provides strong 
evidence for the validity of the model, and our previous deductions that all 13 C02 in our 
experiments was produced via reaction ([2]). 

The results of the model are plotted over the experimental data as open symbols in 
Figure [31 and are summarised in Table [2J At low H irradiation times, the model describes 
well the production of 13 C02, but starts to deviate slightly at longer times. This could be 
attributed to route (jlj) becoming competitive, but in that case, the O3 and O2 experiments 
should deviate with the same trend, while in reality the production of 13 C02 from O3 is 
slightly overestimated while that from O2 is underestimated. Reaction <$5§ is exothermic 



(IMcKinley et al.l 119551 ). so the excess energy of this reaction could be transferred to the 
products, which may then desorb from the surface. The model does not account for this 
possibility, and thus overestimates the 13 C02 yield by the O3 route. 

The model, while a simplistic approach to explaining the surface chemistry, as it ignores 
the complexity of the barrier known to exist in the CO + OH pathway, suggests that the 
relatively low yield of 13 C02 in these experiments results from competition between (j2J) and 
other reactions involving OH. The relation between key reaction rates is: 

kco+H *C kco+OH < kn 2 02+H < kon+H, ko 2 +H, (8) 

indicating that water formation should always dominate the formation of ice species at low 
surface coverages. The overall effective reaction rate of CO + OH was determined to be 
24 times slower than hydrogenation of OH, O3 or O2, and 1.7 times faster than CO + 
H. The relative rate between the hydrogenation of O2 and CO was 40, while between the 
hydrogenation of Q 2 and H2O2, it was 8, consistent with literature values of 31 - 90 and 



3.3, respectively (IMiyauchi et al.l 120081 ) . The factor of two difference between this and the 
previous value of ko 2 +H /kH 2 o 2 +H can be explained by the fact that molecular species such 
as H2O2 were more accessible to H in our (approxima tely 1 ML) experiment s and therefore 



reacted more quickly than in the multilayer regime of IMiyauchi et al.l (120081 ) . 



It is not possible to extract the activation barrier for HOCO formation or its subsequent 
reactions to form CO2 in reaction ([2]) from the model. We can, however, calculate an effective 
barrier to CO2 production by comparison of our modelled relative rates to the activation 
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barrier of 390 K at 12 K for CO + H (IFuchs et al. II2009I ). If we assume first order kinetics and 
a constant pre-exponential factor in both reactions, the effective barrier to (T5]) is 384 ± 40 K, 
where the error is derived from the barrier to hydrogenation of CO. This is the first calculation 
of the effective barrier to (j2J) in the solid phase, and our result is contradictory to that of 
Oba et al (2010) who conclude that, in their study, the reaction proceeds with little or no 
activation barrier, but do not calculate a value. As suggested above, the presence of some 
fraction of excited OH in the beam could produce CO2 in a barrierless reaction with CO. 
Current grain models include grain surf ace activation barrie rs to (T2J) of e.g. 80 K (Y. Aikawa 
2010, private communication) or 176 K (ICuppen et al.ll2009i ). which would seem very low for 
an effective first order rate equation, potentially overproducing C0 2 . All evidence suggests 
that CO + OH is a mo re efficient ro u te to CO2 formation than the non-energetic CO + O 
route, reaction (jTJ, yet iRoser et al. I (120011 ) estimate the barrier to ([TJ is only 290 K, and 



suggest that the reaction proceeds under quiescent clo ud conditions. If this were correct, 
then in comparison to the results of IFuchs et al. I (120091 ) and those presented here, CO + O 
would be more likely to proceed than either CO + H or CO + OH. 

We feel, t herefore, that it is r elevant to address briefly the value of the CO + O barrier 
as derived by IRoser et al. I (120011 ) which, if implemented in gas-grain models, would result 



i n the incorrect pathway to CO2 formation dominating the reaction scheme. IRoser et al 



(120011 ) derived the barrier from an experiment where a water ice cap was deposited on top 
of CO which had previously been exposed to O atoms, assuming that reaction (JTJ occurred 
in the water pores as the surface was heated, and explicitly relying upon the reagents being 
trapped at the surface by the water ice cap . Although other la boratory studies show that 
CO trapping can occur in water ice pores jCollings et alll2003h . observations and models 
confirm that CO freeze-out occu rs after the formation of water ice layers in both molecular 
clouds (IPontoppidan et al.l 120031 ) and protostellar disks (IVisser et al.l 120091 ) . implying that 
a water ice cap is not a realistic mimic of any interstellar ice, including those fo und in 
quiescent regions. In fact, scenarios investigating reaction (jTJ by lRoser et al. I (120011 ) under 
conditions comparable to those prese nt in quiescent mo lecular clouds yielded no CO2, nor did 



subsequent experimental studies by lOba et al. I ( 1201 On . Together with the effective barrier 
to (j2J) of 384 ± 40 K presented here, this suggests the importance of readdressing the value 
of the CO + O barrier implemented in astrochemical models. 



4. Astrophysical Implications 



The reaction of CO + OH is seen to be viable under astrophysical conditions of temper- 
ature and pressure, on silicate and np-H 2 surfaces. Small quantities of CO2 were produced, 
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in competition with other reactions involving OH (for example hydrogenation of OH to form 
H 2 0). Thus, the mechanism CO + OH could be key to explaining the formation of C02( s ) 
at the edges of dark clouds (low Ay), w here C0 2 ( s ) is seen to form concurrently with H 2 
on bare dust grains (IPontoppidan 1 120061 ) . before a large quantity of CO ice is present. 



By modelling the reaction we have determined the empirical relationship kco+H <C 
kco+OH < ku 2 o 2 +H < koH+H, ko 2 +H, where the overall effective rate of CO + OH is deter- 
mined to be 24 times slower than OH + H, and 1.7 times faster than CO + H, indicating 
why H 2 ice is always the most abundant species. 

In dense molecular clouds, gas phase hydrogen is observed to be mainly in the molecular 
form, but atomic H is present at a constant, low abundance (H/H 2 ~ 10 3 



Li fa Goldsmith 



( 120031 )) due to the balance of H 2 formation on grain surfaces and its destructio n by cosmic 



rays. T he abundance of O H increases with density, in line with that of O (see e.g. lHarju et al. 



(I2000[ ); iQuan et all (120081 )). Thus, after the freeze-out of CO, the reaction CO + OH could 
proceed on the ice mantle due to higher abundances of CO on the grain surface. As the 
abundance of OH increases, so does the potential for formation of C0 2 via CO + OH. 
The formation of H 2 via the competitive reaction OH + H will also increase with density, 



and thus C0 2 and H 2 formation in central, more quiescent regions of molecu 



possible. This conclusion agrees well with the postulations of iGoumans et al. I (120081 ) 



ar clouds is 
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Fig. 1. — Chemical network of OH as a schematic diagram. Reactions which occurred in the 
present study are depicted in black, whilst those that didn't are shaded grey. See text for 
details. 
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Fig. 2. — Temperature programmed desorption spectra of Mass 45 ( l3 C02)- The left panel 
shows desorption from a water surface while the right panel, a bare silicate surface. Exper- 
imental data are plotted in their raw form, accompanied by a smoothed version to guide 
the eye. Curves are labelled as follows: solid blue line, experiments A and H (0.5 ML 
O2, ~ 0.5 ML 13 CO, no H irradiation); dashed black line, experiments C and I (0.5 ML 
O2, ~ 0.5 ML 13 CO, 20 minutes H irradiation); dot-dashed red line, experiments G and 
K (~ 1 ML 3 , ~ 0.5 ML 13 CO, 20 minutes H irradiation). Production of 13 C0 2 is seen 
for both starting molecules (O2 and O3) on both surfaces but O3 yields significantly more 
13 C0 2 than O2. There is no discernible surface dependence of the reaction under current 
experimental conditions. 
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Fig. 3. — The evolution of 13 C02 with H irradiation time, presented in monolayers of CO2, 
for all experiments. The data is labelled as follows: red, closed squares O3 + 13 CO on 
np-H 2 0; red, closed triangles O3 + 13 CO on silicate; black, closed squares O2 + 13 CO on 
np-H 2 0; black, closed triangles 2 + 13 CO on silicate. Overplotted (as corresponding open 
shapes) are the results of the kinetic model developed to describe the formation of 13 C02- 
See text for details. 
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Table 1. Experiments performed in this work. 



Experiment 


Substrate 


N(O a ) 


N(0 3 ) 


N( 13 CO) 


t(H) 


(label) 




(ML) 


(ML) 


(ML) 


(minutes) 


A 


H 2 


0.5 




0.6t 





B 


H2O 


0.5 




0.6+ 


10 


C 


H2O 


0.5 




0.6+ 


20 


D 


H2O 




1.6 


0.5 





E 


H2O 




1.4 


0.5 


2 


F 


H2O 




1.6 


0.5 


10 


G 


H 2 




1.1 


0.6 


20 


H 


Silicate 


0.45 




0.45 





I 


Silicate 


0.45 




0.45t 


20 


J 


Silicate 




1.5 


0.13 





K 


Silicate 




1.3 


0.45 


20 



Note. — Species were deposited on the surface in order from left to right, 
apart from those marked f, where 13 CO was deposited first. All species except 
13 CO were deposited using the same beam (see j2]for details.) 
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Table 2. Modelled relative rate constants. 



Reaction 


ki/ko 2 +H 


O3 + H 


1 


O2 + H 


1 


HO2 + H 


1 


H2O2 + H 


0.125 1 " 


CO + H 


0.025 


OH + H 


1 


CO + OH 


0.042+ 



Note. — fRate was a free 
parameter in the model. All 
other rates were fixed, based 
on published empirical val- 
ues detailed in the text. 



